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While two levels systems (TLSs) are ubiqitous in solid state systems, microscopic understanding 
of their nature remains an outstanding problem. Conflicting phenomenological models are used to 
describe TLSs in seemingly similar materials when probed with different experimental techniques. 
Specifically, bulk measurements in amorphous solids have been interpreted using the model of a 
tunneling atom or group of atoms, whereas TLSs observed in the insulating barriers of Josephson 
junction qubits have been understood in terms of tunneling of individual electrons. Motivated by 
recent experiments studying TLSs in Josephson junctions, especially the effects of elastic strain on 
TLS properties, we analyze interaction of the electronic TLS with phonons. We demonstrate that 
strong polaronic effects lead to dramatic changes in TLS properties. Our model gives a quantita- 
tive understanding of the TLS relaxation and dephasing as probed in Josephson junction qubits, 
while providing an alternative interpretation of bulk experiments. We demonstrate that a model of 
polaron dressed electronic TLS leads to estimates for the density and distribution of parameters of 
TLSs consistent with bulk experiments in amorphous solids. This model explains such surprising 
observations of recent experiments as the existence of minima in the energy of some TLSs as a func- 
tion of strain and makes concrete predictions for the character of TLS dephasing near such minima. 
We argue that better understanding of the microscopic nature of TLSs can be used to improve 
properties of quantum devices, from an enhancement of relaxation time of TLSs, to creating new 
types of strongly interacting optomechanical systems. 



I. INTRODUCTION 

At low temperatures many physical properties of amor- 
phous solids are dominated by an ensemble of Two-Level- 
Systems (TLSs/M Measurements of the thermal, elas- 
tic, and dielectric properties of these materials^ have 
been successfully interpreted within the phenomenologi- 
cal model of atoms (or groups of atoms) tunneling in ran- 
dom double well potentials^]. While these experiments 
did not provide a direct evidence for the nature of the 
tunneling objects, the small energies of TLSs were taken 
as corroboration of tunneling atoms rather than elec- 
trons. To contribute appreciably to thermal properties at 
a temperature of hundreds of mK, a significant number of 
TLSs should have energies of the order of GHz, which was 
considered inconsistent with the typical eV energy scale 
for electrons in solids (see e.g. WKB estimate of tunnel- 
ing matrix elements 3 ). However subsequent experiments 
on the effect of TLSs provided evidence that tunneling 
objects carry electric charge and spin, thus casting doubt 
on the picture of tunneling atoms. 

It is now known that TLSs in amorphous solids are very 
sensitive to even weak electric and magnetic fieldsSHS. 
Moreover, TLSs in the insulating layers of Josephson 
Junctions (JJ) were found to be rather strongly coupled 
to JJ qubits. As the energy of the JJ qubit is tuned, it 
comes into resonance with TLSs, leading to anti-crossing 
gaps of up to 30 MHz 9 '^. This strong coupling corre- 
sponds well with the interaction between typical electric 
fields inside JJs and a TLS having a dipole moment cor- 
responding to an e lectro n tunneling through a distance 
of a few Angstrom d 11 * 12 l The aim of this paper is to re- 
solve the seeming discrepancy between the earlier bulk 



measurements, suggesting tunneling of heavy particles, 
and subsequent experiments, indicating the tunneling of 
electrons^HE]. We show that the strong interaction be- 
tween tunneling electrons and phonons leads to signif- 
icant dressing of electrons. We demonstrate that this 
polaronic picture of TLSs not only describes both quali- 
tatively and quantitatively the properties of TLSs as ob- 
served in JJ Qubits, but also provides a natural explana- 
tion for many parameters built into the phenomenological 
model of TLSs given by Anderson et al™ 

Our theoretical model is motivated by recent exper- 
iments with TLSs in JJ qubits by Grabovskij et aZP^, 
which observed dramatic changes in the TLS energy un- 
der the application of mechanical strain to the sample. 
The dimensionless parameter characterizing this change 
K = dlog2E / de, where E is the TLS half-energy and e is 
the strain (which in the simplest case can be understood 
as e = da/ a, where Sa is the change of the lattice con- 
stant a), was found to be consistently in the range of 10 5 . 
Intriguingly, this variation of energy with applied strain 
is in complete accordance with the model of TLS-phonon 
interaction used to describe bulk experiments 

Htls = E z a z + Aa x + j z ea z . (1) 

e here is the strain at the position of the TLS. For TLS 
energies accessible via J J devices (~ 10 GHz), K ~ 10 5 
implies 7 Z ~ 1 eV, which is precisely the value of ^ z 
one deduces from, e.g., bulk measurements of thermal 
conductivity. Furthermore, in the range of strain applied, 
even a complete hyperbolic variation of the half energy 
E = \JE1 + A 2 through its minimum (E z — 0) has been 
observed for some TLSs. 

While such large TLS-phonon coupling motivates the 
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analysis of polaronic effects, one would expect that it 
naturally gives rise to a large rate of dephasmg. This 
stands at odds with another important findingSSEH which 
is that for most, if not all TLSs investigated at energies 
of ~ 10 GHz, the decoherence time T 2 is of the order 
of the relaxation time T\, reaching nearly twice of T\ 
for a few TLSs. As we will show, the bulk of such TLS- 
phonon coupling only leads to polaronic dressing without 
resulting in any anomalously large dephasing. 

The paper is organized as follows. To begin with, in 
Section [TTJ we study the problem of a tunneling system 
in which the principle object is an electron. Typical bar- 
rier height (~ eV) and tunneling length (~ A) result in 
TLS energies at the scale of leV, clearly incompatible 
with observed TLS energies (~ GHz). We show that 
interaction of this confined electron with optical phonon 
modes strongly alters TLS asymmetry E z and suppresses 
A to result in energies in the GHz regime. We then ex- 
plain how coupling to acoustic phonons can result due to 
the modification of the single-electron potential, as well 
as polaronic effects, whereby non-linearity in the optical 
phonon modes provide for a resultant interaction of the 
TLS with acoustic phonons. Our model allows us to show 
why 7 Z ~ 1 eV, explain why j x (coupling of strain e to a x 
neglected in Eq. Q) is insignificantly small, and derive 
a transparent estimate of the spatial density of TLSs and 
the distribution of parameters E z and A. 

Because of strong rcnormalization due to electron- 
phonon coupling, TLS dynamics must be derived fully 
incorporating its polaronic nature. In Section |III[ we 
treat the strong TLS-phonon interaction via a non- 
perturbative variational approach (the partial polaron 
transformation 18 ), that results in an effective Hamilto- 
nian of the TLS that is only weakly correlated with the 
bath of phonons. This then allows us to treat the dy- 
namics of the TLS in the phonon bath under the Born- 
Markov approximation, and analyze the effect of phonons 
on relaxation and dephasing of the TLS. We find, as men- 
tioned, that pure dephasing due to phonons is negligible, 
and that the dominant decay process is via single phonon 
emission leading to an effective Fermi golden rule like de- 
cay rate, which is in quantitative agreement with exper- 
iments. 

From an experimental point of view, an intriguing sig- 
nature of phonon induced decay is that a broad maximum 
(with possible variations due to TLS-TLS interactions) of 
relaxation time T\ is achieved as TLS energy is tuned to 
a minimum by applying strain. This goes against the ex- 
pectation that a fully symmetric double well (achieved at 
minimum E) should be more susceptible to decay from 
a noise that breaks this symmetry (~f z ea z ), but as we 
show, this is really a characteristic of the superohmic 
noise spectrum of strain noised Phonon induced decay 
can also be effectively controlled. In Section |IIID| we 
present an experimental scheme based on band structure 
calculations that show how sizable phononic band gaps 
can be engineered in Al based superconducting qubits by 
appropriate modifications of device geometry, and that 



could result in significantly enhanced lifetimes of TLSs. 

In Section [IV} we discuss the effect of TLS-TLS in- 
teractions on the relaxation and dephasing properties of 
TLSs. We find that other TLSs contribute insignificantly 
to TLS relaxation, but can be expected to dephase the 
TLS in such a way that the decay of the off-diagonal 
TLS density matrix element is gaussian in time rather 
than simply exponential. This altered form of decay is 
a consequence of pure dephasing due to thermally acti- 
vated TLSs that are typically slower than the TLS under 
examination (in the event that the TLS energy is a few 
times greater than the temperature of the system, which 
is usually the case in experiments). Furthermore, at tem- 
perature T ~ 100 mK, the time scale of this dephasing 
process is of the same order as the decay time due to 
phonons of TLSs with energy 2E ~ 10 GHz. This ex- 
plains why most observed TLSs at such energies tend to 
have T 2 of the same order as T\ . 

Finally, we point out that an automatic result of the 
polaronic suppression of tunneling, is the suppression of 
all TLS coupling to the environment through the tun- 
neling (cr x ) operator. As a result, all TLS coupling to 
the environment is primarily through its dipole or a z 
operator, which in the TLS eigenbasis, transforms to 
(E Z /E)(T Z + (A/E)a x This implies that the pure de- 
phasing rate itself can be tuned proportionally to E z /E, 
turning to zero at the TLS symmetry point. At this 
point, TLS decoherence is purely driven by relaxation 
(T 2 ~ 2T 1; and as experiments finoP^}), and our phononic 
crystal scheme becomes extremely viable in significantly 
enhancing TLS coherence properties. 



II. MICROSCOPIC MODEL 




x 



FIG. 1: An electron in a double well potential 

A valence electron typically finds itself sitting in a po- 
tential well structure provided by the electrostatic in- 
teraction with the nearest atoms. Given a disordered 
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spatial configuration of such atoms, it may happen that 
more than one such potential well becomes accessible to 
an electron. In other words, if two such potential wells 
happen to be near enough (say a few A apart), an elec- 
tron can coherently tunnel between such potential min- 
ima, and it becomes important to understand the prop- 
erties of the electron in such a double well potential. The 
double well itself can be characterized by the depths of 
the wells E\ and £2, energies of zero point motion Hui 
and hu>2 in the respective wells, and an interwell separa- 
tion Ax. Given that these structures are created by the 
atomic scale Coulomb interactions, we expect the scale 
of the energies involved to be on the order of an eV. This 
is in contrast to the energy scales used to describe the 
tunneling of atoms, 10 -2 e\S, that is of an order below 
the Debye frequency. For a sufficiently well aligned dou- 
ble well potential, one can work in the tight binding ap- 
proximation of the two lowest lying energy eigenstates in 
either well, neglecting all higher energy levels and arrive 
at the following two level Hamiltonian for the electron in 
the double well structure, 



H = 



- r E 2 + fit* 



(2) 



Here luq is of the order of wi,W2 and T = ^2mV/h 2 Ax, 
m being the mass of the tunneling object, and V is the 
barrier height^. For an electron, with V ~ few eV, and 
Ax of a few A, T ~ 1. In terms of the asymmetry E z = 
Ei + hwi — E2 — Tujj-i and the tunneling strength Ao = 
hu e^ r , one can effectively describe the system by the 
two level system (TLS) Hamiltonian, 



H = E°cr z + A a x 



(3) 



Here E® and Ao are both of the order of an eV. The inter- 
action of the qubit with the TLS is believed to be primar- 
ily an electric field-dipole coupling that can be described 
by Sa z T Xl where a and r are the Pauli matrices corre- 
sponding to the TLS and the qubit respectively. Here, 
S = eE q (Ax) cos</> is the interaction strength that just 
represents the coupling of an electric dipole of charge e, 
length Ax, and orientation angle (the TLS), with the 
qubits electric field whose most significant matrix ele- 
ment is E q , that between the ground and excited states 
(since the electric field of the qubit is proportional to the 
time-derivative of the phase difference across the Joseph- 
son Junction, it is very small in both eigenstates of the 
phase qubit, which are approximately the harmonic os- 
cillator eigenstates in the phase variable). Such estimate 
of interaction was made irP^ and is in reasonable agree- 
ment with the observed gaps in the qubit spectrum due 
to TLSs. Hence, the Hamiltonian that describes the com- 
bined system of qubit and TLS is 



A. Interaction with Longitudinal Optical Phonons 

Before we perform the analysis of dressing of elec- 
trons by high frequency phonons more rigorously using 
a variational approach in Section |III[ to gain some intu- 
ition we consider a simplified model of the interaction of 
the TLS with local Longitudinal Optical (LO) phonons 
whose equilibrium position is strongly affected by the 
state of the TLS. One can write down the Hamiltonian 
as follows, 

1 p 2 

H + H op = E° z a z + A Q <j x + -kx 2 + ^- + Xx{a z +c), (5) 

2 2m 

where one can think of E z as coming from the potential 
created by atoms frozen into an equilibrium position x — 
and H op describes the energetics of the local phonon 
due to interaction with other atoms in the material and 
the electron in the double well. 

Since we expect the energy of the renormalized TLS to 
be in the region of 10 GHz (and as the renormalization 
procedure will yield) which is much smaller than u>d, the 
frequency of the phonons, we expect the phonon, to al- 
ways remain in its ground state, even as the TLS changes 
its internal state. This key assumption will be justified 
later in Section |III| Integrating out the phonon then re- 
quires a minimization of the energy, which is achieved by 
the phonon instantaneously positioning itself according 
to which well the electron is occupying. Since the cou- 
pling of the phonon to the TLS is such that it has different 
energies for different electron positions (due to generally 
nonzero asymmetry constant c), the energy asymmetry 
of the resulting Hamiltonian is also modified by a fac- 
tor proportional to the polaronic shift E p . In particular, 
(x) = —\((T Z + c)/k and E p — A 2 /k. This polaronic shift 
has been estimated to be of the order of an e\^2l in the 
context of Si-MOSFETs, with the material - amorphous 
Si0 2 . 

One must also remember that this procedure also in- 
volves a renormalization of the tunneling strength Ao 
which arises due to modification of T by an additional 
factor of the overlap of the TLS-state-dependent phonon 
wavefunctions. For harmonic oscillator wavefunctions, 
this factor is simply e ~ ( - Xl ~ x ^ t il = e~ Ev ^ D where I is 
the zero-point motion length associated with an oscilla- 
tor of frequency lud- With that, the effective low-energy 
Hamiltonian becomes 



H TT a = (E°-cE„ 



l TLS 

E,, 



A e 



— E p /ujd , 



X 2 /k. 



(6) 



H + Htls-q = E z a z + A a x + So z t x 



(4) 



Note that even though E z and E p both are of the or- 
der of an eV, for some TLSs these can nearly cancel each 
other. Further, since E p /ujd ~ 10, A = e _10 Ao, and 
Ao can as well be of the order of 10~ 5 eV. Thus, renor- 
malization due to optical phonons can in principle result 
in TLSs with energies of the order of 10 -5 eV, which is 
the energy scale at which JJ qubits operate and observe 
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TLSs. Again, the idea here is that the renormalized en- 
ergies of the TLS (determined by many different phonon 
modes collectively) are much smaller than the frequency 
of the phonon (uJd) which means that the phonon al- 
ways remains in its ground state. We want to alert the 
reader that our approach is in contrast with the usual 
idea of the Born Oppenheimer approximation where it 
is the electron that instantaneously assumes the ground 
state wavefunction. 



B. Density and Distribution of TLSs 

As seen above, due to the coupling to LO phonons, 
the asymmetry energy is shifted by cE p . We expect c to 
be a uniformly distributed number of order 1 for optical 
phonons. Therefore, we can obtain a simple estimate of 
the density of the TLSs which have an energy asymme- 
try of about 10 -5 eV. The cancellation of two randomly 
distributed quantities cE p and E®, both of which are of 
the order of 1 cV, to a value of about 10 -5 eV happens 
with a probability of 10 -5 . The total energy of TLS is a 
combination of the energy asymmetry and the tunneling. 
However, due to the polaronic renormalization of the tun- 
neling matrix element (for the typical electron-phonon 
coupling strength in polar insulators, E p /ui£> ~ 10), its 
contribution to total energy is expected to be smaller or 
comparable to the asymmetry contribution's and hence 
does not affect our estimate for the probability. Given 
that it is E z that primarily dictates the energy of the 
TLS, the purely probabilistic cancellation of cE p and E®, 
implies that the energies of the TLSs are uniformly dis- 
tributed in the range of energies E z < E z yielding a 
constant density of states for the TLSs, which is essen- 
tial for obtaining the the linear temperature dependence 
of the specific heat in bulk amorphous solids, and is in 
agreement with Anderson et aiP and Phillips 4 . It is also 
important to emphasize here that such TLSs are expected 
to exist for a wide range of energies (up to eV) , but only 
a small fraction of such TLSs are observable in JJ qubit 
experiments due to the limited range in which the qubit 's 
energy can be tuned. 

What is novel in our approach however is that we can 
go a step further and give a quantitative prediction of 
the density of states of these TLSs. If we take such a 
double-well configurations to occur once in every region 
of about (1(M) 3 , this gives us a density of about 10~ 8 
TLSs per A 3 for TLSs of energy of up to 1CT 5 eV. Using 
the fact the energy of the TLSs are uniformly distributed, 
density of states per unit volume of TLSs obtained is 
10 J _1 m -3 which is indeed in the range of experimental 
estimates (10 45 - 10 46 J" 1 ™" 3 ) of density of TLSs in bulk 
amorphous solidsM 



C. Coupling of TLS to acoustic phonons 

Acoustic phonons can couple to the a z term of the 
TLS by modulating either E p or E® and both forms of 
coupling may be equally important. First, we consider 
the coupling to E z . We can expect that a strain field 
e would modify the energy of either well by a factor of 
E z a t which would give the required coupling strength to 
acoustic phonons of the form ^ z ea z , with j z ~ E* ~ 1 
eV as long as there is an asymmetry between the wells 
to begin with. 

The effect of strain on E p is captured by the Gruneizen 
parameter 77 = d\oguj/de ~ 1 which measures the loga- 
rithmic change in energy of an individual phonon mode 
upon application of uniform strain. Since E p is de- 
termined by optical phonons, which effectively see the 
strain due to long wavelength acoustic phonons as uni- 
form strain, the coupling to phonon fields can also be 
written as E p i]e where, as usual, e is the local strain field 
of acoustic phonons on the TLS. A more detailed deriva- 
tion of this result is given in Appendix [A] 

As we said before, we require cE p and E z to cancel 
each other to give E z in the region of 10~ 5 cV (~ 10 
GHz) to result in TLSs of energies that are accessible in 
JJ Qubit experiments. On the other hand, changes in 
both cE p and E® due to strain are of the order of 1, i.c, 
SE Z /E z e ~ 1 and SE p /E p e ~ 1. Since the changes in the 
two quantities are generically different, the change in E z , 
5E Z — SE Z — c8E p is of the order of eV times strain, i.e. 
5E Z ~ (eV)e. We see that the dimensionless coupling 
parameter K = 8E/(Ee) ~ 5E Z /E z e is essentially a di- 
rect comparison of these two energy scales - eVs and 10 
GHz, which so happens to be the same order of 10 5 that 
is observed in experiments^. We also expect that E z 
and cE p nearly cancel each other for some value of strain 
as it is slowly increased/decreased giving us a minimum 
value of E = ^/A 2 + E 2 Z at some strain value, as is also 
seen in experiments^ (see Fig. [2]). Finally, it is also easy 
to see that the change of A = A e~ E p^ D due to strain 
is much less significant as dlogA/de ~ t]E p /ujd ~ 10, 
which is much less than 10 5 figure that comes from the 
change in E z as a result of the change in E p or E z . This 
explains the absence of a "i x £&z term in Eq. ([I]) 

Before concluding this discussion, we would like to con- 
trast this analysis of the TLS-phonon coupling that arises 
for the electronic TLSs to the case of TLSs resulting from 
the motion of atoms or group of atoms. To find reason- 
able values of the tunneling strength, such as A ~ 10~ 5 
eV, WKB estimates^ require that the atoms live in a po- 
tential landscape with valleys of depths of the order of 
10~ 2 eV. For deeper valleys/barriers, A is quickly dimin- 
ishes to irrelevantly small values. Since phonons interact 
with the TLS by modifying the shape of these valleys, one 
should expect the TLS-phonon interaction of the atomic 
TLS to be of the order of 10 -2 eV times the strain. As we 
know, experiments clearly suggest a much larger value of 
this interaction, it is J z ea z , with ^ z ~ eV. This discrep- 
ancy in the magnitudes of TLS-phonon interaction and 
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the scale of the potential landscape that the constituent 
atoms occupy in the atomic model of TLSs should be a 
cause of concern as both scales correspond to the same 
physics. In our approach, this issue is naturally resolved 
as a confined electron is already associated with eV en- 
ergy scales, and consequently j z ~ eV is well justified. 



III. DECOHERENCE AND RELAXATION OF 
THE TLS BY PHONONS 

TLS decoherence (Ta) and relaxation (Ti) times have 
been measured using the coherent transfer of state infor- 
mation between the TLS and the QubilpH Experiments 
find TLSs with decoherence times T2 almost always of the 
order of T\. Furthermore, as the TLS energy is tuned 
by applying strain, experiments find that at the point 
where the TLS energy is minimum, pure dephasing is 
negligible and T 2 « 2T-P*. As mentioned before, it is 
necessary to take into account polaronic effects in order 
to properly discuss decoherence properties of TLSs. The 
crucial point is that given the large TLS-phonon coupling 
strengths, as found from experiments on both bulk amor- 
phous solids and in JJ Qubit experiments, the phonons 
cannot be relegated to being a perturbative part of the 
environment but actually dress the electron significantly. 
It is this polaronic TLS with its residual interaction with 
phonons that in fact, is nearly immune to pure dephasing 
by phonons, allowing for the observation of TLSs with Ti 
of the order of T\ . Furthermore, we will see that a natural 
result of the renormalization procedure is that the TLS 
can couple to the environment only through its dipole 
operator (<r z in the charge basis). In the TLS eigenba- 
sis, this coupling is of the form - (E z /E)a z + (A/E)a x , 
and thus, as the TLS is tuned to its symmetry point 
E z = 0, pure dephasing becomes negligible and one finds 
T 2 ~ 1T\. We will also show that a properly applied 
Fermi's Golden Rule using the strength of TLS-phonon 
coupling discussed earlier, can correctly capture the mag- 
nitude of the relaxation rate at low temperatures. An 
intriguing signature of phononic decay is that one should 
find an increase of T\ as TLS energy is tuned to a min- 
imum by virtue of the TLS-phonon coupling having a 
superohmic spectrum. 



A. Polaron Transformation 

To treat the problem of dressing of TLS by phonons 
more rigorously than we presented in Section |II A[ we 
use the partial polaron tr ansfor mation approach devel- 
oped by Silbey and Harri^UIMI. advantage of this 
approach is that it does not presuppose any separation 
of time scales between phonons and the TLS, but allows 
for a systematic renormalization of TLS parameters and 
TLS-phonon coupling. The resultant Hamiltonian can 
then be used as an effecti ve mo del to study the dynamical 
properties of the TLS (cf J 23 * 24 * where such a procedure is 



shown to compare well with more rigorous methods) . We 
will also see that this approach indeed consistent with the 
analysis presented in Section [II A| To begin with, we con- 
sider the following Hamiltonian - which involves a generic 
coupling to phonons (acoustic or optical) - 



H = E°a z +A cr x +^2 (°z+Ck) 9k (ak+a'_ k )+^2 u k a\a k . 

k k 

(7) 

This is the Hamiltonian of a TLS coupled to a bath of 
phonons. Note that E z and Ao are still the unrenormal- 
ized TLS parameters, in the eV range as described before. 
Also, hermiticity of the Hamiltonian requires g k = g*_ k 
and c k — c*_ k . Guided by the discussion in the previous 
section, and the form of the coupling deduced from ex- 
periments, the form of the coupling parameter g k is given 
by the term - "f z e(f)a z where e(r) is the local strain field, 
which contains all wavevectors up to kp in the Fourier 
spaced Also, in principle, the strain field has a ten- 
sor structure, and correspondingly so does g k . For the 
purposes of this discussion however, what is important is 
only the power of k in the TLS-phonon coupling, and so 
we do not introduce the full tensorial structure here. For 
more discussion on this, see Appendix [P] 

Now, g k (a k + a'_ k ) = -y z e k giving g k = 
^z.ksjhj (2MNuj k ) where e k is the component of the 
strain field with wavevector k, and MN is essentially 
the total mass of the atoms in the system. Also, c k is 
significant only for high frequency phonons near the De- 
bye frequency (acoustic or optical) which have a wave- 
length comparable to the size of the TLS, that is acoustic 
phonons near the Debye frequency and optical phonons. 
We now transform the Hamiltonian using a generalized 
version of the polaron transformation, 



U = exp 



{a k 



1 



-k) — ( a kVz 
W k 




£7* a k U = a k - 



4 Pt 
—a z 



LJ k 



u k 



f/t al k U = al k ^a z 



(8) 



As is clear from the form of Eq. ([8]), the transforma- 
tion is a matter of simply shifting the phonons accord- 
ing to which state the TLS occupies. The Hamiltonian 
at this step is given in Eq. (Bl|. In what follows, we 



find the optimal values for a k and (3 k by minimizing the 
Bogoluibov-Pcirles upper bound on the free energy (refer 
to Appendix [B] for more details). The result of this is the 
effective Hamiltonian 



9l 



H = (E°-J22c k ^)a z 
^ u k 



Aa x + ^(3fc - Q!fc)(o_fe + ak)<J z 



A (e- 



-2a 



k 

2a 



(e~ 2a ) T )a + + A Q (e 2a -(e 2a ) T )a_, 



(9) 



with the following pertinent equations for a k , E z and A 



G 



= ^ a k /bj k (a k - a_ fc ), 
2E 



«fc = 9k 



1 + — I (AAErtanh/3£coth/3wfc/2 



A = A (e 2a ) o = A exp ( -2 £ 4 coth ^ — ^ 




Wfc. 



Here, E is the half of renormalized TLS energy, E — 
^/E^ + A 2 , A is the "mean field" renormalized tunnel- 
ing strength, while the final two terms in Eq. |9]) serve 
as perturbations. One can see that the factor renormal- 
izing E z has the same form as in Eq. Q, and so we 
suggestively denote factor renormalizing the asymmetry 
by cE p as before. An important result of these equations 
is that renormalization of Ao is primarily determined by 
high frequency phonons. To see this, one needs to note 



that for ui k ^> E, a k 



9k 



which goes as Vk. The 



in the denominator in the exponent of the expression for 
A in Eq. ( 10 1 is naturally cancelled by the density of 



states for phonons, which itself goes as u> k . Thus, due 
to the form of a k , we see that this exponent increases 
linearly in phonon frequency w k . As a consequence, the 
renormalization of Ao is primarily determined by high 
frequency phonons. One can then pull out a factor of 
uj k in Eq (10) for A, replacing it with l)d, and find that 
A ~ A e~^/ WD where E p — ^2 k 2g\j<jj k (compare with 
Eq. (|6|), which is the result we argued for earlier. 

Another consequence of this fact is that since the ther- 
mal occupation of high frequency phonons is nearly un- 
affected by the change in temperature; even though a k 
and A are really part of a self consistency equation, A 
essentially does not change with temperature. This non- 
dependence on temperature is obvious for E z , and this 
result is consistent with the experimental findings^ which 
report no significant change of TLS energies with temper- 
ature. Further, one can evaluate the value of (e 2a ), the 
factor that renormalizes Ao using the coupling strength 
to acoustic phonons as estimated from experiments, i.e., 
j z ~ 0.5eV (and parameters for Al such as density 
p = 2700 kg/m 3 , Debye frequency uj d = 4 x 10 13 rad/s , 
and transverse wave velocity v t — 3000 m/s) and find this 
factor to be of the order of e~ 10 ~ 10~ 5 , that we already 
used above. Thus, as mentioned earlier, one can arrive 
at similar estimates of the suppression of the tunneling 
strength Ao due to known estimates of the coupling to 
optical phonons (as in Section II A), or what is known 



of the TLS-phonon coupling from JJ qubit experiments 
and earlier bulk experiments. 

It is also important to note that this huge suppression 
in fact applies to all couplings that the TLS can have to 
the environment through its a x operator. This can be un- 
derstood by observing that any operator coupled to a x in 
the Hamiltonian Q (before polaron transformation) will 



II A 



acquire a factor (e~ 2a ) after the polaron transformation 
as in eq ([9| . This then implies that the only relevant cou- 
pling of the TLS to its environment is through its dipole 
operator, a z . When we transform into the TLS eigenba- 
sis, this coupling appears as (E z /E)a z + (A/E)a x , and 
thus the coupling to the a z operator in the eigenbasis 
can be tuned as (E z /E). In particular, at the symmetry 
point, where E z = 0, this coupling reduces to zero, thus 
making the TLS nearly immune to pure dephasing, which 
explains why experiments^ commonly find that the de- 
coherence time T 2 reaches twice of T\ at this symmetry 
point. 

Finally, it is important to note that the behavior of 
a k with phonon frequency ui k changes at the scale given 
by the typical frequency of fluctuations in the renor- 
malized TLS's internal state, given by A(A/E). For 
uj k 3> A 2 /E, a k rapidly approaches g k , whereas for 
cu k <C A 2 /E, a k approaches zero. Recall from Eq. |8j 
that a k is associated with displacements of phonons of 
frequency u k in accordance with the TLS's internal state. 
This behavior of a k with uj k thus corresponds to the 
fact that phonons whose frequencies are greater than 
the renormalized TLS tunneling rates, are able to (near) 
instataneously adjust themselves to the TLS's internal 
state (corresponding to the analysis presented in Section 
whereas phonons with much lower frequencies are 
ole to respond effectively to such changes. Thus, even 
though the polaron transformation does not explicitly 
make presumptions about the frequencies of the bath and 
the TLS, the variational procedure automatically incor- 
porates a time scale, associated with renormalized TLS 
tunneling rates, that captures the essential features of 
the response of both slow and fast bath variables. 



B. Energy relaxation 

As mentioned in Appendix |B| the polaronic TLS state 
is a strongly correlated state of the unrenormalized TLS 
and phonon bath, in which the phonon displacements 
are correlated with the TLS's internal state. The op- 
timal values of a k and j3 k describe the mean shifts of 
phonons in correlation with the TLS, and the residual 
interaction with the bath, represents fluctuations beyond 
this mean field shift that has vanishing vaccum expecta- 
tion values. The polaronic TLS, while encapsulating all 
the strong correlations of the unrenormalized TLS and 
phonon bath, is itself only weakly correlated with the 
bath. Thus, one can now apply the Born-Markov ap- 
proximation (see Appendix [C] for more details) to derive 
a master equation describing the decoherence properties 
of the polaronic TLS and attain results for the relaxation 
and dephasing times of the TLS. 

In particular, for relaxation, one finds that even though 
ostensibly Eq. ^ suggests a different rate of decay due 
to change of g k to g k — a k , following the Born-Markov 
procedure results in a (single phonon) decay rate that is 
exactly the same as that obtained from the Fermi Golden 
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FIG. 2: Schematic showing the change in energy E and 
relaxation time T\ of the TLS as uniform strain e 
applied on the system is varied. The Energy E is 
minimal (E = E m i n ) for strain e — e a . At this value of 
the strain, T± concurrently reaches its maximum value 



Rule approach applied to the Hamiltonian with unrenor- 
malized phonon coupling (but with renormalized values 
of E z and A) prior to the polaron transformation. Fur- 
ther, we find that all multi-phonon processes in which the 
TLS decays by releasing more than one phonon, are com- 
pletely insignificant in the temperature range of experi- 
mental interest, and do not affect the decay rate. It must 
also be mentioned that TLS-TLS interactions, which will 



be studied in more detail in Section IV statistically (as 
very few are expected to be resonant with a given TLS) 
contribute insignificantly in comparison to direct phonon 
mediated decay, but when a single TLS's energy is tuned, 
resonances may occur and TLS induced decay may be- 
come significant. The effective relaxation rate of the 
excited state to the ground state can then be found (as 
known from Fermi's Golden rule^) 



(2Ef 
ph 4 



sin 



coth 



E 



(11) 



where sin 9 — A/E and a represents the various phonons 
modes that couple to the TLS. For a TLS of energy 8 
GHz, 7 Z ~ 1 eV, this implies a decay time Ti of about 135 
ns / sin 2 9. For the data on the TLSs presented irP, with 
Ti times of about 400 ns, implies that sin 6* = A/E ~ 0.5. 
It should be noted that it is, in fact, necessary for this fac- 
tor of sin 9 to be significant because the same factor pre- 
cisely enters into the Qubit-TLS coupling strength upon 
diagonalization in Eq. Q and is thus crucial for the TLS 
to produce a visible gap in the qubit spectrum. 

We would now like to point out an intriguing aspect of 
this phonon induced decay rate. Since A remains roughly 
constant as the strain is varied (see section II C I, one can 



glean from Eq. (Ill that the decay rate is actually di- 
rectly proportional to E (recall sin 2 9 = (A/E) 2 ). This 
implies that as the energy of the system is tuned us- 
ing the strain to its minimum value, the lifetime of the 
TLS actually increases to its maximum value (see Fig. 
[2]). This result is unusual because in principle our ex- 
pectation would be for a symmetric double well to be 
more susceptible to decay due to noise that tampers with 
this symmetry than an already asymmetric double well 
(cfP3). In mathematical terms, the matrix element that 
accompanies the decay rate is (A/E) which clearly im- 
plies that decay rate must increase as the energy E is 
lowered. What reverses this trend is the fact that the 
phonon density of states also goes as E 2 besides the fac- 
tor of E that comes from matrix element of the strain 
field. Thus, an increase of Ti as the TLS is tuned to its 
energy minimum is really due to the super-ohmic nature 
of the strain field responsible for TLS decay. This can be 
a simple but effective signature for experiments to check 
for phonon induced decay of TLSs (see also the discus- 
sion in Section IIVI of the effect of TLS-TLS interaction 
on this result.) 

While the above analysis holds true for a single TLS 
whose energy is being modulated by applying strain, a 
similar comment can be made about the distribution of 
the TLSs Ti times vs. energy E. When calculating such 
a distribution, one should be careful to interpret Eq. ( 11 ) 
as really- 
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In principle, we should integrate over the distribution 
of A to find the true distribution of TLSs Ti vs. E. As we 
mentioned earlier, we expect A to be usually smaller or 
at best comparable to the asymmetry energy E z . Under 
this assumption, A does not affect E significantly, and 
we thus expect the scaling of the typical relaxation time 
Ti of TLSs at a given energy E to go as l/E, although 
we expect strong variations due to random values of A. 



C. Decoherence and dephasing 

Pure dephasing due to phonons can be found from the 
Born-Markov approach and the result is mentioned in Eq. 
( |C6[ ). At its largest, the dephasing rate is found to be 
of the order of 1 KHz which is still completely negligible 
in comparison to the decay rate. The pure dephasing 
rate can however also be estimated in a more direct and 
appealing fashion - by estimating the effect of the extra 
phase difference <fi(t) between the eigenstates of the TLS 
generated from the interaction with the phonons - 



4>(t) = 2 
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FIG. 3: (a) The proposed design of the qubit. Note that the holes are drilled in a square lattice pattern at regular 
intervals of 0.24 /jm and the qubit must be fabricated on a substrate no greater than the same length, (b) The band 
structure for k z — 0. T, X, and M represent the points {k x a,k v a) = (0,0), (0,7r),(7r,7r)respectively. The band gap 

extends from cua/vt = 3.57 to ura/vt = 4.23. 



Here 7; ifc e fc = {E z /E)(g k -a k )(a k +a}_ k ), and ji ;k has a k 
dependence due to the presence of a k (the subscript I is to 
indicate the extra factor of (E z /E) in comparison to j z ). 
Using the assumption that the strain field has a gaussian 
distribution implies that < e #(t) >= e"<* (t)2>/2 . Fur- 
ther, assuming that the bath of phonons is uncorrelated 
with the TLS and in thermal equillibrium, we arrive at 
the result 



m 2 ) = 



Hid 



COth {2k^f) {1 - COSUjt) - 
(14) 

We find that the time dependent part has an extremely 
small amplitude of the order of 10 -4 , and performs un- 
derdamped oscillations of the form sin(a>t)/t. Further, 
the time independent part is also small (of the same or- 
der as the time dependent part). Therefore, the phonons 
indeed do not cause any significant pure dephasing, which 
allows for the common observation of TLSs with a deco- 
herence time T 2 ~ 2T\ in experiments^. It is also worth- 



while to note that if we ignore a k in Eq. (14), the time 



independent part of this integral is the essentially the 
same as that in the exponent of the factor that renor- 
malized Ao in Eq. ( 10 ) (setting a k ~ g k for high fre- 
quency phonons). Thus, the polaron transformation acts 
by effectively shifting the effect of constant depahsing 
by high frequency phonons to that of renormalizing the 
tunneling strength. We also already know the estimate 
of this integral to be of the order of 10, which implies 
that prior to this process of systematically addressing the 
role of high frequency phonons in TLS renormalization, 
we would have falsely concluded that the TLS decoheres 
extremely (anomalously such that near instantly, the off- 
diagonal matrix element gets reduced by a factor e -10 ) 



rapidly. 

Finally, we would like to point out that these estimates 
of decoherence are due to interaction with phonons only. 
We do not always expect TLSs to have T 2 ~ 2T\ as 
there will be, in general, other sources of decoherence 
(such as the TLS-TLS interactions analyzed in Section 
IV) coupled to the TLS, making T 2 shorter. 



D. Phononic band gap engineering 

Since TLS decoherence is strongly influenced (with T 2 
being of the order of Xi) by relaxation processes, one may 
expect to get a significant increase in coherence times if 
we are able to increase the lifetime of the TLS excited 
state. It has been demonstrated that due to their fa- 
vorable decoherence and relaxation properties, the TLSs 
can be used themselves as accidental qubits^. There- 
fore, getting longer coherence times could enable new 
types of quantum-computing related experiments. Since 
we find that phonon-induced decay largely explains the 
relaxation rate, the most obvious way of achieving longer 
lifetimes for TLSs is to open up a gap in the phonon spec- 
trum by drilling holes at regular intervals into the Al su- 
perconducting device. Such phononic crystal structures 
have in fact, been recently experimentally realized^. We 
find that a reasonable gap of about 2 GHz can be opened 
up near the 8 GHz frequency marlP^I (Fig. 2b illustrates 
the band gap in terms of dimensionless parameters), by 
drilling holes at intervals of a = 0.24/im in a square lat- 
tice shape, and with a filling fraction / = 0.60. This gap 
however, exists primarily for phonons being radiated in 
the x-y plane (see Fig. 2a). Thus, TLSs can still radiate 
in many other directions even with the phononic crystal 
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structure in place. To make the TLS truly immune to 
decay via phonons, we propose that the fabrication of 
the entire device be performed on a thin film no greater 
than 0.24/im in thickness. This procedure renders the 
phonon of the wavelengths of interest (those that cause 
TLS decay) to become effectively 2D; they now reside in 
a complete band gap, and are in principle, not susceptible 
to decay via phonons. 

As we mentioned earlier, pure dephasing should go to 
negligibly small values (as observed in experiments^) at 
the symmetry point (E z = 0) for TLSs, as their energy 
is tuned by applying strain. At this point, decoherence 
from all noise sources, is almost entirely relaxation driven 
(of which we believe phonons to be the major source 
of energy relaxation), and such an experimental scheme 
may prove to be useful in dramatically enhancing TLS 
coherence times. 



IV. TLS-TLS INTERACTIONS AND 
DECOHERENCE 

Since we find TLS pure dephasing due to phonons to 
be negligible, the role of TLS-TLS interactions in de- 
phasing and perhaps even relaxation must be addressed. 
TLS-TLS interactions have been studied in previous 
works^M!. i n particular, Black and HalperirpHI found 
that by treating the problem classically, neglecting the 
dynamics of the TLSs, one can arrive at an expected 
phonon strairP^ mediated dipole-dipole like interaction 
between the TLSs of the form (in the TLS eigenbasis) - 



Vr = 



Ei Ei 



erf 



5« 



(15) 

where i and j are indices of interacting TLSs, and C%j (f,j ) 
is a material-specific constant that depends only on the 
direction of the vector connecting TLSs i and j. To pro- 
vide for an estimate, Black et aly^* perform rms angular 
averaging over the direction f^, and find for typical val- 
ues of TLS-phonon interaction, phonon velocities etc., 
C„ns ~ 1-6 x 10 -42 Jem 3 . We will use this result of 
Black et a/P51 as the basis for discussing decoherence of 
TLSs, due to TLS-TLS interactions. For completeness 
of this paper, and to establish a connection between the 
polaron approach and previous works, in Appendix ([Pj), 
we show how this interaction follows from the polaron 
transformation procedure in the appropriate limit. 

For the contribution of other TLSs to the decay rate 
of a given TLS, we need to estimate the average in- 
teraction V[ es of our chosen TLS with near resonant 
TLSs. The energy window in which we consider TLSs 
to be resonant is given by max (l/T 1 ,Vf es ), where Ti is 
the relaxation time of our chosen TLS due to phonons. 
We take for example, a TLS of energy 2E ~ 10 GHz, 
and T\ ~ 1/is as found earlier. Using the density of 
states v ~ 10 45 J _1 m -3 , and assuming 1/T X > V[ es , 
the spatial density of TLSs in resonance is given by 



Pres — v(l/Ti), and consequently the average interaction 
strength Vf es ~ C rms p res = C rms v(l/T-C) ~ 0.5 KHz. 
Thus our assumption 1/Ti ^> VJ es is justified. The de- 
cay rate then, is given by (Vf es ) 2 Sx(u = 2Ei) where Sx 
is the correlator of the bath of TLSs that couples to the 
af term of our TLS i in Eq. (15). For large frequen- 
cies, Sx{<*>) is primarily given by the ([txj(oj), crj(— + ) 
correlators^ of the bath TLSs, that are Lorentzians cen- 
tered at u) = 2Ej with width I/T2J, where T2J is the de- 
coherence time of TLS j. For the case of resonant TLSs, 



Ej ~ Ei, the decay rate is given by (V[ es ) 



2,/ 



10- 



KHz which is significantly smaller than the contribu- 
tion to decay by phonons. Here we used the exper- 
imental observation that for most TLSs in the energy 
range of 10 GHz, as our example TLS, T^j is of the or- 
der Tij ~ 1/MHz. Thus, statistically, we expect that 
TLS decay due to other TLSs is negligible. Neverthe- 
less, we would like to point out that when a selected 
TLS's energy is tuned by applying strain, it is bound 
to come into resonance with some TLSs. Coupling be- 
tween such TLSs depends on the distance between them, 
so it should have strong variations from one resonance 
to another. At such TLS resonances we expect to find 
considerable variation in the T± in the range of strain, 
where the difference of TLS energies is comparable to 
their coupling. These rapid resonance variations of the 
relaxation rate can make it difficult to verify the smooth 
trend T\ oc 1/E profile that we obtained earlier. 

We now look at the contribution of other TLSs in the 
dephasing of our TLS. We expect only TLSs that are 
thermally active (of half energy Ej below the tempera- 
ture T) to contribute to dephasing. Assuming that the 
temperature of the system is lower than the energy of 
the TLS we are examining (under typical experimental 
conditions T ~ 100 mK is a few times smaller than TLS 
energy 2E ~ 10 GHz), these thermally activated TLSs 
are expected to be slower than our TLS (recall from Eq. 
( 12 ), decay times typically decrease with decreasing TLS 



energies). The dephasing due to such slow TLSs can be 
obtained only after averaging over many experimental 
runs, each of which would correspond to a different but 
near static configuration of the bath. The flipping of such 
slow TLSs thus, merely shifts the energy of our TLS, at 
the scale of the typical interaction energy yp° n - res _ As 
we discuss below, dephasing due to such TLSs is qualita- 
tively different from a bath that operates at time scales 
faster than our TLS. 

For the general case of a gaussian bath, one finds the 
pure dephasing induced decay of the off-diagonal density 
matrix element of the TLS is given by^ 



MO) 



:exp(- 2 



1 f du> . ,sin 2 (u;£/2) 
2n Sz{UJ > ( W /2) 2 



(16) 



where Sz{u) is the symmetrized bath correlator in the 
Fourier domain, S z (w) = 2 f^dte^ {[V z (t),V z (0)] + ). 
Here, V z is sum total of the interaction terms from all 
TLSs j 7^ i that couple to the erf operator of TLS i in 
Eq. (15). At low frequencies, the major contribution to 
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Sz(oj) comes from the correlators J^j ([ cr |( w )i <J j{~ UJ )] + ) 
that are given 

' (17) 



uPTlj (cosh Ej/kbTf 

Thus, only TLSs with half energy Ej < k^T contribute 
to dephasing, which is why we consider only thermally 
activated TLSs. Given the 1/r 3 nature of the interac- 
tion between TLSs, nearest neighbor TLSs are expected 
to most strongly affect TLS decoherence properties, and 
make a dominant contribution to the net bath correla- 
tor Sz{u). Motivated by these observations, to bring 
out the difference between pure dephasing due to a slow 
bath and a fast one, we consider a bath with correlator 



Sz(u>) = A 2 r/(r 2 



), which is a telegraph noise with 



a switching rate L and interaction strength A, and we 
evaluate the exponent in Eq. ( 16 ) in the two regimes, 



1/r > t, appropriate for a slow bath and 1/L -C t for a 
fast bath, where the time f that we're interested in, is set 
by the time scale at which the TLS under examination 
decoheres (say around 1/is). Now, the scale at which the 



integral decays in the exponent of Eq. ( 16 ) is set by ei- 
ther Sz(w) which decays at the scale V or (sm 2 (u;t/2))/w 2 
that decays at the scale 1/f. For a short-correlated bath, 
we find (cr_(t)) ~ exp (— (A 2 /T)t) . For a slow bath, with 
T <C 1/t, we find instead (<r_(t)) ~ exp (— A 2 t 2 ) wherein 
the only time scale is given by I /A. In the present case, 
we are in the regime where the TLS bath's switching rate, 
set by T - 1/Tij < 1/t ~ lfis. That is, the TLS bath 
is slower than the TLS under examination. 

Thus, in the case of dephasing of a given TLS by other 
TLSs, the only scale that decides the additional contribu- 
tion to dephasing is given A, or in this case, yp on ~ res an d 
the decay is gaussian, that is, exp (-(y r non_rea ) 2 t 2 ) . To 
estimate yp on - res i we fi rs t note that the density of ther- 
mally activated TLSs is given by p n on-res ~ vT, where 
T is the temperature and, consequently yp on ~ res ^ 
C rms vT ~ 0.5 MHz, at T = 100mA'. Since this is of 
the order of linewidth of our TLS, we expect strong vari- 
ations in the values of I/T2 whereby the nearest neigh- 
bor thermally active TLS strongly influences the deco- 
herence property of a given TLS. When this nearest 
neighbor distance is shorter than the typical distance 
r typ ~ (1/i/T) 1 / 3 - 100 nm for T = 100mA, we ex- 
pect interactions with TLSs will dominate the dephasing 
of our TLS. However, if this distance is larger than rtyp, 
we expect TLS decoherence to be set by phonons, with 
T 2 ~ 2Ti. Nevertheless, from the above analysis, it seems 
reasonable to expect, that for most TLSs of energies of 
the order of 10 GHz, T2 will be at least of the same order 
as Ti, which is what one finds in experiments. 

To complete this discussion, we point out that the dif- 
ferent nature of decoherence due to phonons (via decay) 
and TLS-TLS interactions (via pure dephasing) allows us 
to simply add up their corresponding contributions to the 
decay of the off-diagonal term of the TLS density matrix. 
In a single experimental run, the TLS bath shifts the fre- 
quency of our TLS slightly, simply modifying oscillation 



frequencies in a Ramsey experiment. The gaussian form 
of decay due to TLS-TLS interactions appears upon av- 
eraging over many experimental runs. Hence, we suggest 
that to analyze such experimental data, one should care- 
fully fit the decay with a more general 2 parameter fit 
of exp {—A 2 t 2 — Bt) , where B should be approximately 
l/(2Ti). It is also important to note that the interaction 
A = v™ n ~ res cx [E*/E)i for the TLS % under observa- 
tion. This means that at the point of TLS energy min- 
ima (Ef = 0, when tuned by applying strain), one should 
achieve minimal dephasing. At this point, decoherence 
is truly driven by phononic relaxation, that can be min- 
imized via the fabrication scheme described in Section 

eho 

Again, we remind our reader that this tuning of pure 
dephasing to zero at the symmetry point of the TLS holds 
very generally in our theory. Since all the coupling of the 
environment to the TLS via its a x operator is strongly 
suppressed (by the same large factor that the tunneling 
is suppressed by), the TLS's coupling to the environment 
is expected to be primarily through its dipole operator, 
represented by a z in the charge basis. Consequently, the 
interaction of the TLS with the environment in its eigen- 
basis must of the form - (E z /E)a z + (A/E)<r x , which 
means that as we approach the symmetry point of the 
TLS, where E z = 0, the pure dephasing rate becomes 
entirely negligible, yielding a decoherence time T2 —> 1T\ 
at this point. As far as we know, experiments^ are in 
complete agreement with this statement. 



V. CONCLUSIONS AND OUTLOOK 

We conclude with the understanding that tunneling 
systems observed in bulk amorphous solids, as pro- 
pounded by Anderson et alW and Phillips^ are actually 
electrons dressed by optical frequency phonons tunnel- 
ing between atoms rather than atoms or group of atoms 
tunneling in double well structures. This reasoning is 
motivated by the fact that coupling strengths suggested 
by experiments both on bulk solids and in JJ Qubits are 
too large to ignore their impact on the renormalization 
of the TLS parameters. Further, once we assume that 
it is indeed an electron that is the tunneling object, ex- 
planation of the coupling to acoustic phonons becomes 
transparent, the density of such TLSs can be estimated 
to the correct order of magnitude as seen in experiments, 
and the coupling of such TLSs to the JJ Qubit is also 
easily explained. 

Due to renormalization of TLSs by optical frequency 
phonons, the decoherence and relaxation properties are 
also significantly altered, and in particular pure dephas- 
ing due to phonons becomes insignificant. TLS dephasing 
is then given primarily by interaction with other TLSs, 
while the low temperature relaxation rate is given pri- 
marily by phonon induced decay. In particular, we find 
that pure dephasing due to the bath of TLSs leads to 
a gaussian in time decay of the off-diagonal density ma- 
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trix element, with an effective time scale given by the 
average interaction of the TLS with other thermally ac- 
tivated TLSs in the bath. For T = 100 mK, we find this 
time scale for pure dephasing to be of the order of the de- 
cay rate (due to phonons) of TLSs of energy in the GHz 
range, explaining quantitatively, why for most TLSs at 
such energies, experiments find T 2 to be of the order of 
Ti if not twice, with T x ~ 

We also provide an explanation for the existence of a 
minima in the energy of the TLS as a function of strain, 
and predict that phonon induced decay must produce the 
intriguing signature of a maximum of relaxation time T\ 
at this symmetry point contrary to usual expectations. 
The renormalization procedure we present also has the 
important consequence that any coupling to the environ- 
ment must be through the dipole (a z ) operator in the 
TLS charge basis. This directly leads to the practically 
important result that the dephasing rate of the TLS must 
become insignificantly small as it is tuned to its symme- 
try point by the application of stain and explains why 
experiments consistently observe that Ti ~ 2T\ at this 
symmetry point^. An important outstanding issue is 
that of explaining the temperature dependence of the 
decay rate of the TLS; we believe this to be a result 
of coupling not only to phonons, but also to the other 
parts of Qubits environment, including non-equilibrium 
quasiparticles 3 ^ in the superconducting Al layer. 

Finally, we presented a schematic of an experimental 
setup that should result in suppressing phonon-induced 
decay of TLSs, enhancing their lifetime, and conse- 
quently decoherence times, significantly. This experi- 
mental scheme becomes more viable as TLSs are tuned 
to their energy minimum as strain is varied, since we ex- 
pect dephasing by other TLSs at this symmetry point to 
be minimal. While we don't expect this to improve the 
coherence properties of superconducting devices them- 
selves, (for who m these TLSs are a major source of 
decoherenc d 34 * 37 ^ ^) ; this scheme could have an impor- 
tant impact on not only improving the viability of TLSs 
for use as quantum memories^, but should also allow us 
to understand better the other sources of noise that de- 
cohere the TLSs, once the influence of phonons is signifi- 
cantly reduced. A better understanding of these aspects 
of the TLS should eventually help improve the practica- 
bility of both superconducting devices and optomechani- 
cal resonator^] where these TLSs significantly limit the 
coherence times and Q-factors respectively, as well as en- 
gineer exciting new nano-mechanical devices with strong 
non-linearities of phonon modes. 
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Appendix A: TLS-acoustic phonon coupling induced 
by phonon nonlinearity 

The origin of Gruneizen parameter lies in non-linear 
phonon coupling, Hnl = r yx 2 e(r). The form of this term 
is based of the following intuition - the first generic non- 
harmonic term has to be cubic in the phonon fields; it 
must involve two optical phonon field terms whose fre- 
quencies nearly cancel each other so as to give a resul- 
tant acoustic frequency, and that the coupling to acous- 
tic phonons must be via a derivative of the acoustic 
phonon field. In effect, Hml serves to modify the stiff- 
ness constant k of the optical phonon mode, which al- 
lows us (using the definition of the Gruneizen parameter 
i] = dloguj/de, k — raw 2 ) to arrive at the following rela- 
tion - 7 = rjk. Thus, we can write the following Hamil- 
tonian 

H + H op + Hml = E z cr z + Ao<j x + 
p 2 1 

+ -k (1 + rie(r)) x 2 + Xx(a z + c). 
2m 2 

(Al) 

Eliminating the optical phonon mode now follows the 
same procedure as before - the only difference being a 
modified value for k, and we finally obtain 

Htls + H ac = E° z <r z + Aa x - cE p / (1 + rje(r)) a z . 

(A2) 

Expanding the above in T)e(r) gives us the desired 
acoustic phonon-TLS coupling, and hence reproduces the 
Hamiltonian in Eq. Q. Since r]e(r) < 1 (at T = lOOmK, 

^\A 2 ( r )^ < 10~ 8 ), all terms beyond the first order in 
e(r) can be safely neglected. 

Appendix B: Polaron Transform 

We begin with the general Hamiltonian coupling a TLS 
and a bath of phonons as in Eq. Q and apply to it the 
polaron transformation U mentioned in Eq. (fSJ) . Then 
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H = WHU is given by 
H = E' C + E' z a z + Aa x + ^ (a k + al k ) d z (g k - a k ) 

k 

]T u k a\a k + A a+ (e~ 2a - (e" 2 ")) 

k 

+ A a_ (e 2a - (e 2a )) + ^(a fe + al k )(c k9k - fa). 

k 

(Bl) 

Here a are the pauli matrices in the [/-transformed basis. 
The values of E' c and E' z are 



El 



E 



o^fcCfcfffc 


2 a* k g k + \<*l\ 




U k 


Wfc UJ k 





E> = El - 



E 

k 
2a ' 



Wfc 



, a fcCfcfffc 
Wfc 



A = A <e 2Q ), 



fc 



(B2) 



Note that the reality of these terms is guaranteed by 
the conditions on the hermicity of the Hamiltonian and 
the unitarity of the polaron transformation (which re- 
quires, a k = a*_ k , and fa — P*_ k ). Now, from the form 
of U in Eq. ([8]), one can see that the polaron transforma- 
tion is simply a unitary transformation that allows us one 
to go from a basis of independent TLS (a) and phonon 
variables, to the new TLS variable (a) that corresponds 
to a correlated motion of phonons and the unrenormal- 
ized TLS. The philosophy then is to somehow optimize 
these correlations such that the free energy of the pola- 
ronic TLS system, which now is only weakly correlated 
to the residual bath, is minimized. To do so, one intro- 
duces the set of parameters a k and fa , that describe the 
amplitude of phonon displacements correlated with the 
TLS's internal state and provide a variational space over 
which to minimize the free energy of the polaronic TLS. 
The optimal values of a k and fa thus, represent the mean 
phonon displacements in correlation with TLS state, and 
the residual coupling with the bath represents fluctua- 
tions about this mean displacement, that has vanishing 
vacuum expectation values. 

Given the polaron-transformed Hamiltonian H in Eq. 
(Bl I, we need to estimate the minimum of the free en- 
ergy. But first, we must complete the squaring o f th e 
phonon field to eliminate the last term in Eq. ( Bl ) . 



This is brought about by the transformation - d k = 
a}_ k — [fa — c k g k )/u) k . In fact, this transformation elimi- 
nates fa (replacing it with c k g k ) from the transformed 
Hamiltonian completely leaving us with a variational 
state dependent only on a k . It is important to note that 
this last transformation leaves the coefficients of cr+ and 
<t_ unchanged in Eq. (Bl ) and these can now be treated 



as perturbations to the phonon and TLS parts of the 
Hamiltonian, as their expectation values in the thermal 



ensemble of phonons is zero. In the process of separat- 
ing these these final terms, we arrived at the mean field 
estimate for the renormalized tunneling, A = A (e 2a ) . 

To perform the approximate minimization of the free 
energy, we minimize the Bogoluibov-Peierls upper bound 
on the free energy (which here is simply the free energy 
of the TLS and phonons decoupled from each other) F = 
E c — (l/(3)lnZ + F p h- Z is partition function of the TLS, 

e -PE' + e fSE' with E i = ^ E ,2 + A ^ and p ph is thc frec 

energy of the phonons which is independent of a k . The 
result of this is the following equation on a k , 

{ A \ 2 IE 

K-Pfe)+ (-E?) — a k twahl3E catia.0u k /2 = O. (B3) 
V El uj k 



Solving this leads to thc main results in Eq. (10). Also, 
note that since g k is real for this case, so is a k and we 
consequently drop the usage of extra symbols such as * 
and II. 



Appendix C: Calculation of decay rate 

The first step in this procedure is to diagonalize the 
TLS part of the Hamiltonian in Eq. ([9| . The Hamilto- 
nian can then be represented by 





H = 


where 


V = 


§E(3* 

fc 


v? = 


WE, 
2\E 




§E(^ 

fc 


X = 


(e- 2a < 



H = Ea z + V?a + + Vycr- + V Q a Zl 



(CI) 



(i ^ v- , A A^ t 



An A 
2 E 



2 E 



1 A X 



1 (E z 
E 



- 1 A Xt 



<e- 2 «» 



(C2) 



The time evolution of the density matrix in thc Born- 
Markov approximation^ is described by the equation - 



dp 



dsTr B [Vi{t), [Vi{t - s),p s (t) ® p B }} 



(C3) 

Expanding the RHS of Eq ( |C3| ) out, and keeping only 
the slowly fluctuating terms (the rotating wave approxi- 
mation), we arrive at the estimates for T\ and Ti - 



1 



1 _ 1 
T~ 2 ~ 27\ 



dse l2Es (v?(s)V 1 (0)\ + h.c. 

Vi(8)V}(0)) + h.c, 
(V o (s)V o (0)) . 



dse 



-i2Es 



4i?e 



(C4) 
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More explicitly, 
l - = 2ir f^) (9k ~ a k f (2n k + 1) 6(u k - 2E) 



■ 8tt (~) AV — (jt- a k ) (2n k + 1) 5{u k - 2E) 

\ E J 

e aEs ds (e G ^ - 1 




+ 1 \Re 



1 Re 



1 Re 



- 1 .Re 



e aEs ds (e- G ^ - 1 
e- aEs ds (e G ^ - l) 



-i2E. 



: ds (e~ a ^ - l) 
(C5) 



1 _ 1 

% ~ m 

where, 



2A J 



Re 



ds (e~ G ^+e G{ ^-2) 
(C6) 



G(s)=4Vf^ ) ((n k + l)e-^ s +n k e^ s ). (C7) 

Adding up the processes involving only 1 phonon oper- 
ator in Eq. (C5| gives the result in Eq. (111. As for the 



pure dephasing rate, T$ — 1/T 2 — l/2Ti being negligible, 
a quick way to see this is to realize that if we expand the 
exponentials in Eq. (C6|), the first non vanishing term 
comes with two Bose distribution functions, proportional 
to n(uj)(n(uj) + l) (with the same frequency) which decays 
rapidly for frequencies at which a k is appreciable. 



Appendix D: TLS-TLS interactions using Polaron 
Tr ansfor mat ion 

To get a full description of the TLS-TLS interaction, 
one needs to include the full tensorial structure for the 
strain field that couples to the TLS. In particular, one 
has e a p{r) = (1/2) {d a up{r) +dpu a (r)), where u a (r) is 
the phonon displacement field at position r. To include 
this full description, we also have to consider the 3 differ- 
ent phonon operators or equivalently the transverse and 
longitudinal modes of the phonon fields, and append g k 
with tensor indices af3. Again, for the purpose of this 
paper, to get an estimate of the interaction strength be- 
tween TLSs, we avoid this detail that does not change the 
power law of the interaction but simply introduces more 
tensorial structure, and show how one can attain the ex- 
pected 1/r 3 interaction via the polaron transformation 
method. 



We now consider the Hamiltonian of two TLSs (identi- 
fied by indices 1 and 2 at positions and f resp.) coupled 
to the phonon strain field as before, 

H= E V°°t + A?o? + X>iK + cJ)(a* + al fc ) 

i=l,2 k 

+ 92^2 + c k 2 )(a k + al k y k - f + J2 u k a{a k . (Dl) 

k k 

In what follows, we will include e lk ' r in g k . 

We now introduce the Polaron transformation, that 
allows for phonons to shift according to both TLS's states 

U = exp ^~X!( a * - a -k)^ a i a i + a 2 a 2 + Pk)j ■ 

(D2) 

For the Unitarity of U, we require a k = (a i k )*, 
[3 k = P*_ k - This is supplemented with the Hermitic- 
ity of the Hamiltonian Eq. (Dl|, that demands that 
Si = (fli )*• The transformed Hamiltonian under the 
polaron transformation is now 

H = E C + E{o\ + Elal + A^af (e" 2 " 1 ) + A°af (e~ 2 " 2 ) 
+ E^lul + J2( V i°l + ^2 2 )K + a_ fc ) 

k 

+ £ Lo k a\a k + Y(a k + a)_ k )X k 
k k 

+ £ A?o-+(e-*" - (e- 2 -)) + *>i(e +2at ~ (e- 2a % 

(D3) 



i=l,2 



where X k = f3 k — c\g\ + c\g\. Again, before we can 
estimate the free energy, we need to eliminate the free 
linear term in the phonon fields, which requires shifting 
the phonon fields by X k , as defined in Eq. (D3|. As 

c 252 ' 



in Section [Bj this amounts to setting f3 k = c\g1 1 



thus eliminating f3 k from the variational procedure. The 
expressions for the terms in our Hamiltonian now reduce 
to (with X k = 0) 



\c\gX + c\g\\ 



UJ k 

k V , ("2 - 25 2 fc ) 



+ K) 



k „ {a\ - 2g k ) 



u k 



UJ k 



E* = E r-J2 2 < 



k W 



u k 



l^coth/?w fe /2 , 
\ /. ^ J 



2{a\)* 



2{a k Y 



k 



LO k 



U! k 



V k = 9 



k a k . 



(D4) 
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Thus, we find again that Ef's are independent of the 
variational parameter oq, and that Aj is determined 
by the phonon displacements due to each TLS indepen- 
dently. However, if Ej ^ 0, then oq's are coupled and 
the renormalization of A< is indeed influenced by other 
TLSs. 

In principle, one must now minimize the free energy 
considering the spectrum of the Hamiltonian Hmf — 
Ec + Y^i^t a i + Ajof -I- Eja{crl, which is not a sim- 
ple task. It is also important to remind oneself that 
we cannot interpret Ej itself as the complete interaction 
strength between TLS 1 and TLS 2. The polaron trans- 
formation approach merely sets up an effective Hamil- 
tonian in which the phonons have been shifted to their 
mean field positions. The remnant interaction with the 
phonon fields have zero expectation values, or are fluctu- 
ations which in principle can give an interaction between 
the TLSs. These in general involve considering the full 
dynamics of TLSs and phonons, and cannot be replaced 
by a simple effective Hamiltonian. Such a task is beyond 
the scope of this paper. 

The results of Black et al. however can be achieved by 
considering the regime in which spins do not have any 
dynamics, i.e, in the limit A — > 0. In such a scenario, the 
spins are essentially static and the system is classical. 
The only variables that depend on a^'s are now Ej and 
E c . Regardless of the temperature of the system, mini- 
mization of free energy requires simply setting oq = gf . 
This can be seen by noting that d ^ky = (<?i — oi\) and 

= (.92 - that is > onl y proportional to a\ - gf. 
This simplification has the important consequence that 
the remnant interaction with phonons are also completely 



killed, — which allows us to understand Ej as the 
complete interaction energy. In this scenario, we find - 



Ej 



{glYg k 2 



*9x 



k 



U) k 



U) k 



(D5) 



Now, for 



g\ 



9k and 

r y z k^H/2pu> k as before, we find that the sum over k 
in Eq. ( D5 ) is simply a fourier transform over some- 
thing that is independent of any power of k, that is, 



g\ 



ik.r 



g k e^- with g k 
find that the sum over 



Er 



Ak.r 



This may seem peculiar in 



that, ordinarily this would imply a <5 function interaction, 
but it must be noted that in 3 dimensions, the fourier 
transformation of any T(f)/r 3 interaction is also inde- 
pendent of any power of k but in general involves depen- 
dence on the direction of k. Suggestively, if we include 
anisotropies of the form (ih\.k) and (ifi2-k) (with h\ and 
hi defining dipole axes for TLSs 1 and 2 respectively) in 
g\ and g\ respectively, we arrive at the classical dipolc- 
dipole interaction with the magnitude ^/(Airpv 2 ^) and 
the appropriate tensor structure. 

For the purpose of this paper, it suffices to show that 
the polaron transformation can be used to yield the effec- 
tive interaction between TLSs in the limit that A — > 0, 
that is, upon neglecting the effect of TLS dynamics on 
phonon displacements. For the more complete result, the 
reader should look at the paper by Black and HalperirPSl, 
where upon obtaining the complete TLS-TLS interaction, 
they perform rms averaging over the angular structure of 
the interaction to arrive at an estimate of C rms in Eq. 



(15). Finally, to arrive at the result in Eq. (15) one sim 



ply rotates the basis to diagonalize the inc 
Hamiltonians in Eq. (D3). 
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